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ClY Abstract 

O l' We study cosmic-ray anomaly observed by PAMELA based on Eq inspired extra U(l) model with 

S4 x Z2 flavor symmetry. In our model, the lightest flavon has very long lifetime of O(10 18 ) second which 
i-C ■ is longer than the age of the universe, but not long enough to explain the PAMELA result ~ 0(1O 26 ) 

sec. Such a situation could be avoidable by considering that the flavon is not the dominant component of 
, dark matters and the dominant one is the lightest neutralino. With appropriate parameter set, density 

^ ■ parameter of dark matter and over-abundance of positron flux in cosmic-ray are realized at the same 

I time. There is interesting correlation between spectrum of positron flux and Vmns- No excess of anti- 

proton in cosmic-ray suggests that sfermions are heavier than 4 TeV and the masses of the light Higgs 
r**j~ ' bosons are degenerated. 
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1 Introduction 



Standard Model (SM) is successful theory of gauge interactions, however Higgs sector is not examined well. 
Therefore mass matrices of leptons and quarks are not well understood. Many unsolved puzzles of SM are 
left in these sectors; that is, e.g., why is the structure of mixing matrix of leptons (Maki-Nakagawa-Sakata 
matrix, Vmns) very different from that of the mixing matrix of quarks (Cabibbo-Kobayashi-Maskawa matrix, 
Vckm), especially why is the mixing angle #23 maximal? Why is neutrino mass far smaller than those of 
other fermions? Why do generations exist? 

We also find a problem in cosmology. In modern cosmology, the existence of the dark matter is clear. 
Recent cosmic-ray observation of PAMELA suggests that the dark matter decays mainly into leptons with 
very long lifetime [1][2][3]. Such a particle is not included in SM. 

Separately from these puzzles, there is hierarchy problem why electroweak scale is much smaller than 
Planck scale. One of the solutions is to introduce supersymmetry (SUSY) [5]. However minimal super- 
symmetric standard model (MSSM) does not satisfy the solution, because we must fine-tune //-parameter in 
superpotential of MSSM, which is much smaller than Planck scale in order to realize appropriate electroweak 
symmetry breaking. This is called //-problem. 

Another problem of MSSM is proton stability. The R-parity forbids baryon number violating trilinear 
terms in superptential, however does not forbid quartic terms like E C U C U C D C , LQQQ. Such interactions 
reduce the lifetime of proton to unacceptable level [6] . Therefore the R-parity does not help the explanation 
of proton stability. The problem of proton lifetime of supcrsymmctric model is one of the most essential 
point in understanding generation structure. 

With the motivation to solve flavor puzzles and hierarchy problem, we introduce new three symmetries. 
At first, we introduce non-Abelian discrete flavor symmetry S4 x Z 2 , in order to explain that the mixing 
angle #23 is maximal [7] [8] [9] [10] [11]. Because Vckm and Vmns are verv different, it is expected that the 
representations of quarks and leptons are also different. Next, we introduce U(l)x gauge symmetry which 
forbids //-term [12]. Then, several new superfields must be introduced due to gauge anomaly cancellation 
condition; those are extra Higgs (H U ,H D ), singlet Higgs S and exotic quarks (g,g c ). The extra Higgs 
bosons couple only to leptons, which induce the difference between Vckm and Vmns- Moreover, the exis- 
tence of exotic quarks is important to understand the meaning of generations. Finally we introduce U(l)z 
gauge symmetry. Due to the anomaly cancellation condition, right-handed neutrino (RHN) superfield N c is 
introduced, then the smallness of neutrino mass is realized by seesaw mechanism. The two new U(l) gauge 
symmetries and standard model gauge symmetry Gsm = SU(3) C x SU(2)w x U(l)y can be embedded in 
E 6 as Gsm x U(1)x x U(1)z C E 6 , then MSSM and new superfields consist 27 of E e representation. With 
appropriate assignment of superfields under the flavor symmetry, the stability of proton is realized, which 
plays the most important role in the flavor symmetry Thus we can understand that the generation structure 
is the new system to stabilize proton [8] . 

The new symmetries which are introduced above may also solve dark matter problems. As three U(l) 
gauge symmetries include R-parity, lightest supersymmetric particle (LSP) is a candidate for dark matter. 
The positron flux observed by PAMELA is produced by the field which induces RHN mass and decays into 
leptons [13]. In this paper, we show our model is consistent with experimental results of dark matter. At 
first, we define our model in section 2. The estimations of relic abundance of dark matter and positron flux 
are given in section 3. Finally, we give conclusion of our analysis in section 4. 



2 S4 x Z2 flavor symmetric extra U(l) model 
2.1 Gauge symmetry 

We extend the gauge symmetry from Gsm to Gi = Gsm x U(1)x x U(1)z C Eq, and add new superfields 
S,g,g c ,N c which are embedded in 27 representation of Eq with quark, lepton superfields Q,U C , D c , L, E c 
and Higgs superfields H u ,H D . In order to break U(l)z gauge symmetry, we introduce Gsm singlet $ and 
$ c . The gauge representations of these superfields are given in Table 1. After the gauge symmetry breaking, 
as the R-parity symmetry 
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Table 1: G*2 assignment of superfields. Where the x, y and z are charges of U(l)x, U(l)y and U(l)z, and 
Y is hypcrcharge. The extra U(l) charges x and 2 are given by x = \q^p + \q x and z = — + \q x , where 

<?</> = & \JIQtP and <7 X = 2^(9^ are charges of two U(l)s in SU(5) x U(l)j, x f/(l) x C 5*0(10) x {7(1)^, c 

remains unbroken, LSP is the candidate for dark matter. The invariant supcrpotcntial under the gauge 
symmetry G2 is given by 

W = Y U H U QU C + Y D H D QD C + Y E H D LE C 

+ \SH U H D + kSgg c 

+ Y N H U LN C + Y M $N C N C 

+ M$$ c $ + Vl QQg + y 2 g c U c D c + y 3 gE c U c + y i9 c LQ + y 5 gD c N c , (2) 

where first line consists of trilinear terms in MSSM. Second line generates effective /x term A (S) H u H D by 
radiative symmetry breaking of U(\)x- Third line generates RHN mass term Y N ($) N C N C by radiative 
symmetry breaking of U(l)z and gives small neutrino mass by seesaw mechanism. Fourth line consists of 
unwanted terms which cause the problems such that the mass term M$$$ c prevents $, $ c from developing 
vacuum expectation values (VEVs) and the trilinear terms of exotic quarks destabilize proton. Note that 
Higgs superfields are extended to three generations. Generally, extra Higgs doublets cause the problem of 
flavor changing neutral currents (FCNCs). 

2.2 Flavor symmetry 
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Table 2: S4 x Z 2 assignment of superfields (Where the index i of the ^4 doublets runs i = 1,2, and the index 
a of the 64 triplets runs a — 1, 2, 3.) 

In ordre to explain maximal mixing angle #23, we introduce S4 x Z 2 flavor symmetry. This symmetry 
solves the problems of supcrpotcntial defined in Eq.(2) at the same time. If we assign g,g c , $ c to S4 -triplets 
and the other superfields to singlets or doublets, then M$, yi,....5 arc eliminated. As a result, supcrpotcntial 
is given by 

W' = Y U H U QU C + Y D H D QD C + Y E H D LE C 
+ XSH U H D + kSgg c 
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+ Y N H U LN C + Y M <PN C N C 

4- _^_<j><j><j> c <j> c 

M P 

+ TJ2^ C QQ9 + ^^ c g c U c D c + ^^ c gE c U c 

1V1 p Ivl p Ivl p 

+ JJT^ C 9 c LQ+^^ c gD c N c + ---, (3) 

where the dots ■ ■ • are higher order terms. As the potential of $ and $ c is lifted by non-renormalizable term 
<f><f><f> c <f> c , $ and $ c have very large VEVs along the D-flat direction of ($) = (<f> c ), where 

From the constraints on the lifetimes of proton and exotic quarks (see appendix C), the condition 

must be satisfied [14]. From this condition, V ~ 10 12 GeV is required. This value is realized by potential 
minimum condition as Eq.(4), when we take a ~ 10~ 2 . The prediction for RHN mass is given by 

M R ~ 10 12 GcV, (6) 

which gives an appropriate neutrino mass. As the VEVs of $ and & c break not only U(l)z but also S4 x Z 2 , 
we call them flavons. 

2.3 Maki-Nakagawa-Sakata matrix Vmns 

The maximal mixing angle #23 of Vmns is realized by the assignments that H u , H D , L, N c are 2 + 1 of S4 
and E c is 1 + 1 + 1' [7]. In order to reduce the number of parameters, we assign S and $ to 2 + 1. If we 
assign Q 7 U c and D c to ^-singlets, then quarks do not couple to S4 doublet Higgs and FCNC is suppressed. 
The flavor representations are given in Table 2. The leading order superpotential is given by 

W SiX z 2 = WL + Wg + Wn + Wg + W^, (7) 

W L = Y 2 N [H? [LiN c 2 + L 2 N{) + HZiLi.Nl - L 2 N^)} 

+ Y 3 N H^L 3 N^ + Y 4 N L 3 (H"NC + H% V 2 C ) 

+ YfEKHfLi + H?L 2 ) + YfEZHgLz + Y 3 E E^{H[ ) L 2 - H?L X ) 

+ ^Y 1 M ^ 3 (N^ + N^) + ^Y S M ^ 3 NINI, (8) 

W Q - Y¥ H 3 QiUj + Y^H 3 QiDj = 1,2,3), (9) 

= XiSsiHYH? +H^H 2 D ) + X 3 S 3 H^H 3 D 

+ XiH^iSxH? + S 2 H. 2 D ) + \ 5 (S 1 H? + 82^)^, (10) 

VF 9 = kS 3 ( gi gi + g 2 g^ + g 3 gl), (11) 

W* = 2^*l[(*f) 2 + (*§) 2 + (*§) a ] 

+ ^ {2V3^-f 2 mr - mf] + ($ 2 - + ($§) 2 - 2^) 2 ]} . (12) 

We give the parameter set to realize the maximal mixing angle of MNS matrix. We define non-negative 
VEVs as 

(H?) = (H") = ±v u , (HV) = v' u , {H? ) = {H?) = ±v d , (H») = v' d , 
(Si) = (S 2 ) = -j=v a , (S 3 ) = v' s , 
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= Vi, ($ 2 ) = Jvf + v\ + v\ - vl ($ 3 > = V, 



<*?> = Vir +w?> &*> = \Jt +v1 <^} = y^+*< 2 « v). (13) 

In Wz,, without loss of generality, we can define Y% 4 , 3^ 3 , to be real and define the phase of Y 3 N as 
Y 3 N = \Y 3 N \e lS . We define mass parameters as 

Mi = Y^V, M 3 = Y 3 M V, 

m%=Y 2 N v u , m» = \Y 3 N \v' u , m\=Y^v u , (14) 
m\ = Yfv d , mi 2 = Y 2 E v' d , m\ = Y 3 E v d . 

Using these parameters, the mass matrices of charged leptons and neutrinos are given by 

n\ —mk \ ( TO, m, 



Mi = I m[ m\ , M D = -m£ 

V2m l 2 / V TO 4 m 4 \/2e l V 

Mi 
Mi 
M 3 




sin 

sin 6^ cos Q v 



where 



(15) 



The following neutrino mass matrix is generated through the seesaw mechanism 

/ p\ p 2 p4 \ 

M v = MoM^Mjj = p\ , (16) 

V p\ + e 2l Vl / 

where 

m 2 m 4 m 3 f-\ rr\ 

p2 = vm' PA = 7m' P3 = vm- (17) 

The mass eigenvalues and diagonalization matrix of charged leptons are given by 

V^MfMfr = diag(mlml,m 2 T ) = ((m l 2 f,(mi) 2 ,(m[) 2 ), (18) 

1 I _1 M 

Vi = — ( 1 1 I , (19) 

and those of neutrinos are given by 

V£M v Vu = diagie^-Wm^e^+Wm^m^), (20) 
/ -sin 6^ e^cosOu \ 
V u = 1 I . (21) 

\ e^cosBu sm6 v 

From Eq.(19) and Eq.(21), we obtain the MNS matrix as follows 



Vmns = V?V v P v = -^=\ sm6 v -e* co&9 v 1 I P„, (22) 



P v = diag(e- lWl -^ )/2 , e - 4 (^+«/ 2 , 1). (23) 



Here it is worth mentioning that the lower bound of (0.04 <)#i3 was shown by the recent experiment reported 
by T2K [15] at 90 % C.L., which could give a severe test to our model near future. 
From the experimental bound [18], we impose the condition 

tan(9, y = — , ml 2 -ml ± =8.0 x 10 _5 (eV 2 ), to 2 2 - m 2 ^ = 2.5 x 10~ 3 (cV 2 ), (24) 
v2 
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on the parameters, then the phase <f> is given by 



r cos ( 



0.361, 



Fixing the VEVs as 

v u = 10, v' u = 155.3, v d = 2.0, 
and the charged lepton masses as [8] [19] 



P2 

Pi 



i 



- 1.75GeV, m 2 = 487kcV, 
Yukawa coupling constants are given by 



Y 



E 



0.875, Y 3 E = 5.15 x 10" 2 , 



77.8 (GeV), 



103MeV, 



Y E = 6.25 x 10" 6 . 



For the RHN mass parameters, we assume 

V 



10 12 GeV, 



Y{ 



M 



Y, 



M 



1. 



(25) 

(26) 
(27) 
(28) 
(29) 



In order to investigate model dependence, we give two sample parameter sets A and B which arc defined as 
follows 



0.361 

: 01 = 2 = 0.0°, S 

= 1.80 x 10 _2 eV 



90.0°, 



p\ = 14.47 x 10~ 2 cV, p\ = 13.78 x 10~ 2 cV, 



5.24 x 10 _2 eV, 



= 5.31 x 10" 2 eV, 



1.80 x 10" 2 eV, 



m 2 = 4.24GeV, m u 3 = 12.0GeV, 



11.7GeV, 



Y 



N 



B 



= 0.424, 

1.000 

68.84°, 



Y 3 N = 0.077, 



41.55° 



Y 4 N = 1.17, 



62 = 41.15°, (5 = 89.805°, 



(30) 



p\ = 5.03 x 10 _2 eV, p 2 = 10.04 x 10 _2 eV, 



p\ = 5.03 x 10~ 2 eV, 



7.08 x 10" 2 eV, 



7.14 x 10 _2 eV, 



5.03 x 10" 2 eV 



m l z = 7.09GcV, 



mX, = 10.02GeV, 



7.09GeV, 



Y 



N 



0.709, Y 



N 



0.065, 



Yf = 0.709. 



(31) 



Generally, multi-Higgs model causes FCNC problems, however our assignments do not cause such problems. 
In the lepton sector, the interactions between charged leptons and Higgs bosons are given by 



Yfr c 



V2 
HP + HP 



v/2 



+ l T 

+ It 



HP + Hj 
V2 



Hi 



Y 2 E H»e% 



D 



HP 



V2 



(32) 



which do not contribute tor— s>e + 7,/i— >e + 7 processes. Because Li has accidental 5*2 symmetry such as 

{HP, HP) -j. {HP, HP), {l^p c )^{-l^-p c ), (33) 

t — > p + 7 process is also not induced. Note that this S2 symmetry is not the symmetry of whole theory, as 
the symmetry is violated in neutrino Yukawa couplings and flavon superpotential . 

In the quark sector, as the quarks couple only to H 3l Hp, Higgs mediated FCNCs are not induced. In 
the basis that quark mass matrices are diagonal, the superpotential is written as 



W Q = YtHg(Q 3 y(UZy + Y c H^{Q 2 )'{Ui)' + Y U H% {Q^' {U<[)' 
+ Y b HP{Q 3 )'{D c 3 )' + Y s H 3 D {Q 2 )'{D c 2 y + Y d H E \Q 1 )' '{D^) 1 \ 

From here, we fix top, bottom and charm masses as [8] [19] 

Y t v' u = 172.5, F 6 ^ = 2.89, Y c v' u = 0.624 (GeV), 



(34) 
(35) 
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and assume exotic quark mass as 

fc<=2000 (GeV), (36) 

in order to forbid the decay of lightest flavon into exotic quark pair. Then we fix the values of Yukawa 
coupling constants as 

Y t = 1.12, Y b = 0.0371, Y c = 0.00405, k = 1.0, v' s = 2000GeV, v s = 200GcV. (37) 



2.4 Higgs sector 

Higgs potential is given as follows, 

V = V F + V D + V A + V m 2, (38) 
Vp = \\1S3H1 + X^S\H 3 \ + I A1S3H® + \c,S2HP^ 



+ I ^3S 3 H 3 D + \^(SiHi + S 2 H. 

+ J \\S3H1 + \iSiH^\ + J \1S3H2 + X4S 2 H 3 \ 

+ \\ 3 S 3 H% + \ 5 (SiH? + S 2 H%)\ 2 
+ + X§HY H 3 \ + \X±H 3 + X5H2 H 3 ^ 

+ ^{H^H^ + H^H^) + X 3 H U 3 H 3 D \ 2 , (39) 
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2 \ 2 a 2' a 1 / 2' 

x ' A 

+ ^(-2|^| 2 -3|^| 2 + 5|5 a | 2 ) 2 , (40) 

V A = -X^SaiH? H? + H$ iff) - \ 3 A 3 S 3 H% 

- XJUH^S^H? + S 2 H?)-\MS 1 Hl J + S 2 H%)H? + /i.e., (41) 
V m , = -^ 3 |^| 2 + ^(|^f| 2 + |^| 2 ) + ^ 3 I^| 2 + ^(l^fl 2 + l^l 2 ) 

- m| 3 |5 3 | 2 +m|(|5i| 2 + |S2| a ) 

[m 2 BU (H^(HY + H?) + m% D (H°yiH? + H?) + m 2 BS (S 3 )\S 1 + S 2 ) + h.c] , (42) 

where we can define Ai,3 ^5 to be real without loss of generality, and we assume all the soft SUSY breaking 
parameters are real to avoid complex VEVs. 

he soft S4 x Z 2 breaking terms; m BU ,rn BD ,m BS , violate accidental 0(2) symmetry of Higgs potential 
and fix the VEV directions (Eq.(13)) to realize 6*23 = 45°. This potential has S 2 symmetry such as 

H% O H%, H? o H?, Si^ S 2 . (43) 

Minimizing this potential, we get mass matrices of Higgs bosons. The results are given in appendix A. In 
the same manner, we add soft S4 x Z 2 breaking terms in flavon sector to avoid domain wall problem [20] . 



3 Dark Matter 

Here we show that our model is consistent with cosmic-ray observation of PAMELA. Decaying dark matter 
scenarios with Non-Abelian discrete flavor symmetries have been done by Ref. [16]. 

3.1 LF decay width 

We assume that the candidate for decaying dark matter is the lightest flavon (LF). In the six flavon super- 
fields; $ , $a(a = 1, 2, 3), only one linear combination is super-heavy and the other five superfields have TeV 
scale masses. As LF cannot decay into other flavons, it has very long lifetime. Due to the non-renormalizablc 

1 We would like to thank Refree for the suggestion. 
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interactions with light particles, LF becomes unstable dark matter. Among the interactions, the source term 
of RHN mass 



(Y N H U L) 2 (Y N H U L) 2 ( $ 3 \ 1 [H U L\ f $ 3 . 
%= ,vM/.x .r =o^ — 1-17 (44) 



rN H u L) 2 (Y N H U _ 
2Y M (V + <Z> 3 )~ 2Mr V" V r y'~2"^V " 7 V" V" 

is the unique interaction to emit leptons without emitting quarks [21], where w is VEV of H u . We estimate 
the positron flux using this interaction. Due to the factor 1/v, the Higgs which develops the smallest VEV 
gives the largest contribution to LF decay. Therefore we can neglect the contribution from , because 
v u <C v' u as one can see from Eq.(26). This effect is impotant to suppress weak boson emission. Due to the 
enhancement factor m^p/vu, LF decay width is dominated by 4-body decay as follows 

T(LF -+v + v) < T(LF — >• i/ + Z + H + ) ~ T(LF v + 1 + W + ) 

< T{LF + l + H + + H + ). (45) 

From the spectrum of positron flux observed by PAMELA, we assume 

m iF =4TeV. (46) 

If we assume all sfermions which couple to LF are heavier than 4TeV, the other interactions do not contribute 
to LF decay. The interactions which contribute to LF decay is given as follows 

C 2v = ±C LF <t> LF {[(HY) (HV) + {H^f{H^f]{v e v e + r 2 ^ + r 2 v T v T ) 

+ 2V2r{v, - v T ){HYf( H V fv e - v/2r(i/„ + Vr )[(H?) (H?) - {H^f{H^f]u e ) , (47) 

Civ = -\c LF 4> LF {2[{H^)\H^) + + {H^f{H^) + ]{eu e + r 2 v^ + r 2 v T r) 
+ 2r 2 [(HVf{HV)+ - (HV)°(H?) + }(v,T v r y) 
+ V2r[(HV)°(HY)+ + {HVf(H¥) + ][{v» "r)e + i/.fc r)} 

V2r[(HYnHY) + -(H^nH^)+}[(^ + u T )e + ^ + r)}} 7 (48) 

C-21 = 1 -C lf Hf{[{H¥) + {HV)+ + (H^)+(H^)+](ee + r 2 W + r 2 rr) 

+ 2V2r( M - T ){H?) + (H?)+e V2r(fi + t)[(H")+ (H?)* - (H") + (H^)+]e} , (49) 

C » = Till' (50) 



where e is the flavon mixing parameter which is defined by 

V2 



*3 = ^^, (51) 



where <f>L F is LF field. 

Using Eq.(47)-(49), the LF decay widths are given as follows 

r a „ = T 2P = (6 + 10r 2 + 12r 4 )r , (52) 

r lv = r & = r e + r |i + r r = (2 + 8r 2 + 8r 4 )r 0> 

r e = r g = (2 + 4r 2 )r , r M = = r r = i> = (2r 2 + 4r 4 )r , 

r 2( = r 2r = r 2e + r 2M + r 2T + r eM + r eT = (8 + i2r 2 + i6r 4 )r , (53) 
r 2e = 8r , r 2A1 = r 2T = 8r 4 r , r eM = r er = 6r 2 r , (54) 

r lepton = r^+^ = (10 + 20r 2 + 24r 4 )r =r anti _ lcpton , 
r total = 2(r 2!/ + r^ + r 2; ) = 2(16 + 30r 2 + 36r 4 )r , (55) 
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where wc classify the final states only by charged lepton flavor e, /i, r. The rates of lepton flavor emitted by 
LF decay are given by 

9 + 8r 2 4r 2 + 10r 4 , . 

Pe = 7T-T^r-^KU> P^=Pr= . , | =t-t . (57) 



9 + 16r : 



20r 4 



9 + 16r 2 + 20r 4 ' 



Anti- lepton flux depends on Majorana phase <\> through Eq.(25), such as e-dominant for r < 1 and (/i,r) 
dominant for r > 1 (see Fig. 1). For each parameter set, LF lifetime is estimated as follows 



A 
B 



total 



total 



3.72 x 10 n e" 2 sec, r~\ • , , = 1.17 x 10 12 e" 2 sec, 
' anti-lcpton ' 

= 2.13 x 10 12 e- 2 sec. 



= 7.01 x 10 n e- 2 sec, T 



anti-lepton 



Hereafter we assume e < 10 to avoid extinction of LF. 




(58) 
(59) 



Figure 1: 



3.2 Relic Abundance of LF 



At early stage of the universe, Savon multiplets are produced through U(l)z gauge interaction [13]. Since 
we assume that reheating temperature is low enough to avoid gravitino over-production as Trh < 10 GeV 
[22] , this interaction is never thermal equiliburium. Therefore we assume non-thermal production of flavons 
and boundary condition tilf(Trh) = 0. 

For the chiral multiplets (4>l, U(l)z gauge interaction is given by 



C 



i )z = ig z A^ * [^Pi,L^tpi,L + *i^*t 



(60) 



from which we calculate production cross sections of flavon multiplets (tp, $). From Eq.(13), the U (1) z gauge 
boson mass is nearly equal to 16g z V. 

As all produced flavon multiplets decay into LF finally, LF number density is given by 



(61) 



where 5 is the number of light flavon superfields, are number density of one flavon multiplet and 

Nlf ~ 0(1) is LF production rate which means how many LFs are produced per one degree of freedom of 
flavon multiplets. The Boltzmann equation for Ulf is given by 

2480N LF CT 8 , 



3Hn F p 
C 



21 



1 



2,2V 2 



from which we get 



n LF h 2 



m LF s h 2 



(2tt) 5 

15 x 2480 x 21m P N LF ( 1 



(62) 
(63) 



2tt 2 (341.25) x 30.67(2tt) 5 yi2V 2 



rp3 

1 RH 



5.06 x 10~ 9 N LF 



( T RH 
\ 10 5 GeV 



(64) 
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where [18] 





T 2 


H 


= 1.66-y/g* , 




mp 


.</» 


= 341.25, 


mp 


= 1.22 x 10 19 GeV, 


so 


= 2890/cm 3 , 


Pc 


= 1.05 x 10 4 /i 2 cV/cm 3 



For TpH < 10 7 GeV, LF does not dominate dark matter (Jl^ph 2 <C £1dm^ 2 = 0.11), thus other dark 
matter should be considered as we will discuss later. Such multi-component dark matter is discussed in 
[26]. Although the number density is very low, the short lifetime of LF enables us to explain cosmic-ray 
observation. The effective lifetime of LF is defined as 

r eff = r anti-lcpton ' (66) 

and the following values are obtained for each parameter set 

4.2 x 1(T 7 , (67) 

6.6 x 1(T 7 . (68) 
Eq.(67) and Eq.(68) are satisfied for example, if we put N LF ~ 1, T RH ~ 10 5 GcV, e - 1CT 3 . 



A 



Trh \ 



r cS = 6.0 x 10-sec, e<N LF (^^^ j 
B : r cff = 7.0xl0 25 sec, e 2 N LF (^g^) 



The positron flux from the decay of LF is calculated as 

*(E e+ ) = ^f^— [ dE'G e +(E e +,E') V Pi^r, (69) 

47T TTlLFToS J n f—', , dE' 

where v e + is the velocity of the positron, G e + is the Green's fuction which is expressed in [3] , pt is expressed in 
Eq.(57) and dNg e + j dE' is the fragmentation function produced from the decay of I to e + . The fragmentation 
function is calculated by using the event generator pythia [27] and the result is shown in Fig. 2. We can 
evalutate the positron flux from the decay of LF by using the fragmentation function. The results for each 
parameter set A and B are shown in Fig. 3. 




Figure 2: The fragmentation function calculated by pythia for parameter set A (left) and B (right). 



From the gamma-ray observations [4] , the constraint for r-flux is given by 
, « r lepton 10 + 20r 2 + 24r 4 in26 

( T -)cff = —£T— T C ff = 8r 2 +20r 4 T cff > 2- 1 X 10 SeC ' ( 70 ) 
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which is estimated for each parameter set as follows 

A : (r T ) cff = 5.6 x 10 26 sec, 

B : (r r ) cff = 1.4 x 10 26 sec. (71) 

As the parameter set B is severe to satisfy Eq.(70), small r model is favored. This is the new information 
about ncurtino sector extarcted by cosmic-ray observations. 

3.3 Higgs decay width 

No excess of anti-proton flux in cosmic-ray constrains the species of the particles emitted by LF decay [23]. 
As the weak boson Z, and the chargino decay mainly into quarks, LF should not decay into these 
particles so much. The weak boson emission is suppressed by factor (v u /niLF) 2 and the chargino emission 
channel is kinematically closed for heavy sfermion scenario. In order to forbid the weak boson and the 
chargino emission from Higgs boson decay, we assume light Higgs scenario. 

In the Higgs potential Eqs.(38)-(42) and mass terms of the neutralinos and the charginos 

C D -iV2{H^[g 2 X^ + 2,g y X Y ~2g x X x ]h u a 

- iV2(H°y{g 2 \$Tf - 3g y X Y - 3g x \xK - iV2(S a ^[5g x X x }s a 

- \*h>4*2 - \m y Xy\y - \m x \x\x + (W H )g2 + h.c. (A = 1, 2, 3), (72) 



we assume the parameters as follows 



32 = 0.652, g y = igy (g Y = 0.357), g x = -^=g Y = 0.073, 

Ai = 0.065, A 3 = 0.4, A 4 = 0.398, A 5 = 0.75, 
M x = M 2 = 200, M Y = 180 (GcV), 
A x = 0.0, A 3 = Ai = A 5 = 1.0 (TeV), 



2 _ n O/ioQ ™2 
BU 



0.0408, m% D = m* BS = 0.02 (TcV^). (73) 



Mass matrices of neutralinos and charginos are given in appendix A and the values of mass eigenvalues and 
mixing matrices are given in appendix B. 

We consider only the mass eigenstates which dominate H^ 2 such as 

01(91.50), 04(121-96), ^(152.48), ^(112.12), ^(145.53), 

(if 1 ± )'(90.70), (i? 3 ± )'(130.02) (GcV), (74) 

where <f>[, p[, (H^)' are ^-odd and the others are even. As S 2 forbids interaction (f>[ZZ and 4>[ is not 
emitted through Z* — > Z + <$>\ , LEP bound run > 114.4GeV is not imposed on (/)[. As the masses of these 
Higgs bosons are well degenerated, they do not emit weak bosons or charginos. 
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The ncutralinos into which these Higgs bosons can decay are two singlino dominant neutralinos 

V [ (41.92), 774(44.55) ( GcV ), (75) 



where r( x is SVodd and LSP. S^-even neutralino 774 can decay into r( x through 774 — > rj[ + fx + f without 
emitting quarks. As SVodd Higgs boson can not decay into quarks, we consider only the decay of S^-even 
Higgs bosons. 

The decay widths of 4 , 5 , p' e due to the Yukawa interactions 



C D Y c (HV)°cc c + Y b (H»)%b c + YfKH^e^ + {Hffefa] + h.c. 
are given as follows 





— > c - 


he) 


= 3.63 x 10" 


t 2d 


r(0 4 


-> 6- 


h&) 


= 1.54x10" 


t 2d 


r(0 4 - 


-> 


- f) 


= 1.10x10" 


t 2d 




-> c - 


he) 


= 1.61 x 10" 


-5-p 
1 2d 


r(0' 5 


-> 6- 


h&) 


= 2.68 x 10" 


" 4 r2d 


r(0' 5 - 


-^t4 


- f) 


= 4.90 x 10" 


t 2d 




— > c - 


he) 


= 1.91 x 10" 


" 7 r2d 




-> 6- 


h&) 


= 4.96 x 10" 


t 2d 


r(p 6 - 


-^t4 


- f) 


= 6.89 x 10" 


-8 r 
1 2d 



(76) 



(77) 
(78) 
(79) 
(80) 
(81) 
(82) 
(83) 
(84) 
(85) 



where T 2 d is 2-body decay width of scalar. The interactions with the ncutralinos rj^,^, 



C D 



-X, {S 3 [(h«)°(h? )° + + s 3 [(H??(h? )° + (i^) ^) ] 

S3 [(^) (i/f) + (^)°(i/ 2 D ) ]} 

A3^3(^)°(^) + -a(2*3W) + -sC^W )°] 

A 4 {(HV)°[ Sl (h? )° + S2 (^)°] + (>#) [Si(>if )° + &(#)°] 

(^)°[ Sl (iff)° + S2 (iff) ]} 

A 5 + *a(^)°] + (^)°[Si(^)° + S 2 {h u 2 f] 

(h°)°[si(HY)° + s 2 (HVf]} 



i 
i 

iV2^T Sj[5g x X x ]si + h.c 



2 92X 



-g Y X Y - 2g x X x 



^9yXy — 3g x Xx 



{hi)° 
(h?)° 



-(0.062004 + 0.12705 + 0.0299ipe)»?4'74 
(-0.071604 + O.1170' 5 + 0.0142^)77^; + h.c. 



give 



r(0 4 -+ v + v ) 
r(0 5 ^^ + 77) 
r(Pe -> 77 + 77) 

which dominate the decay widths of 4 , 5 , p' 6 . 

The decay widths of (H^)' due to Yukawa interactions 



208 x 10~ 4 r 2 d, 
864 x 10- 4 r 2d , 
3179 x 10~ 6 r 



2d- 



C D -Y c {H")+sc c -Yf[(Hf ) )-e c v 1 + (H?)-e'iis 2 }+h.c. 



(86) 

(87) 
(88) 
(89) 

(90) 
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are given by 

T((H-)' -> s + c) = 1.95 x 10- 7 r M , (91) 
T{{H-)' y t -\- v T ) = 2.48 x lCT 6 r 2d . (92) 

From these estimations, (H^ 1 )' decay gives dominant contribution to anti-proton flux. 

As one charged lepton emission from LF decay accommodates one charged Higgs emission at even rate 
3 

-l 

^eff = 27 - 4 ^eff- ( 93 ) 



of (H^Y and (H^ 1 )' , the quark flux is estimated as 



r quark 

For each parameter sets, we get 







0.195 



2.48 + 0.195 



A : r quark = 1.6 x W sec, (94) 
B : T quark = 1.9 x 10 27 sec, (95) 
from which the spectrum of anti-proton flux is given in Fig. 4. There is no inconsistency in anti-proton flux. 



x=1 0x10^ sec 




0.1 1 10 100 

T [GeV](m DM = 4.0 TeV) 



Figure 4: 



3.4 Relic Abundance of LSP 

Finally we estimate the relic abundance of LSP. The interactions between r][ , rj'^ and Z are given by 
C D ^i7 A1 (-5 Al -iG 1 Z Al 7 5 )Vi + ^47 M (-^-»G4^75)V'4 



where 



+ iGih^f^Z^j + iGif^frZ^f, (96) 



d = 0.00823, G 4 = 0.0119, (97) 

0.172, G(e R ) = f^M= = 0.200, (98) 



G(v L ) = - v;y V = -°- 372 > (") 



9y 




V9y + 


si 




9l _ 


2 










9l 


9 Y 




3V9y 


±92 



3.114, G{u R ) = -^=M= = -0.257, (100) 
6V3y + 92 

-0.057, G(d R ) = -4jM= = 0.315. (101) 

6 V9y + 92 
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These interactions give dominant contribution to annihilation of 77^ , 77^ . As Gi <C G(z^l), the contributions 
Z — > rjrj to Z-decay width is negligible. Therefore LEP bound m > 46GcV is not imposed on Tj[,r]'^. The 
relic abundance of the neutralino is calculated by the formula 



xf = In 



O.Q955mprrii(a + 66/xf) 



Qh 2 



(g*x F ) 

8.76 x 10- 11 



(* = 1,4), 



g* xf 



(a + '3b/x F )GeV 2 



(102) 
(103) 



where mi = 41.92GeV, m 4 = 44.55GeV 



2c/ 

7T 



^§( Gtfi ,- G (/ H ), 



1 - 



(104) 



4 mf 



■■in 



m t (G t /2) 
4to 2 — ml- 



[G 2 (h) + G 2 (f R )] 4- 



2m / 



m 
m 



(105) 



c/ is color factor such as c/ = 1 for SU C (3) singlet, c/ = 3 for triplet, and mz = 91.2 GeV. For the 
approximation mf/mi = 0, these coefficients are given by 



016,1 + Ctr,l 
Q.&.4 + a-r,4 



a 1 
"4 

61 = 1.53 x 10~ s GcV 

64 = 4.23 x 10~ 7 GcV~ 2 



5.01 x 10~ n GcV~ 2 , 



1.22 x 10~ 9 GcV~ 2 , 



(106) 
(107) 
(108) 
(109) 



The relic abundance of r][ is given by 



and that of 774 is given by 



9* 

n x h 2 



9* 

Xf 

n A h 2 



75.75, 

22.32, 
1.88 GeV, 
0.106, 



72.25, 
25.52, 
1.75 GeV, 
0.0052. 



(110) 



As ?7 4 is converted into 77^ , relic abundance of LSP is given by 



(flGDAlh 2 ) = flih 2 

which realizes density parameter of dark matter. 



n 4 h 2 = 0.111, 



(111) 



(112) 



4 Conclusion 

We have considered dark matter based on S4 x Z2 flavor symmetric extra U(l) model. The results are as 
follows. There exists appropriate parameter set to realize relic abundance of dark matter and positron flux 
observed by PAMELA at the same time. The dominant component of dark matter is LSP and the origin 
of positron flux is given by the decay of LF which generates the mass of RHN. There is deep connection 
between PAMELA observation and neutrino mass. The long life time of LF results in large RHN mass and 
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the spectrum of positron flux depend on Majorana phase <f> in Vmns- Therefore, cosmic-ray observation 
gives new information about the structure of neutrino mass matrix. 

From the fact that there is no excess of anti-proton flux in cosmic-ray, we can guess about the particle 
spectrum. As sfermions can decay into quarks, weak boson and charginos, LF must not decay into those 
particles, which sugests sfermions are heavier than 4TeV. Although this is also favorable from the viewpoint 
of the FCNC constraints, experimental verification becomes difficult. However experimantal verification of 
our scenario is not imposssible. From the fact that Higgs does not decay into weak boson or charginos, we 
can expect that Higgs boson is light and degenerated, therefore the examination of the mass spectrum of 
Higgs boson is possible. 



A Mass matrices 

Neutral CP even Higgs boson 

r U <fiu,a 



H, 



V2 



-C D 



1 



VP. a 

V2 



D% (a = 1,2,3), 



V2 




(i,3 



1,2, 



,9), 



M l,3 



ML 



AI 



4.4 



M 



M 



4,5 
2 

4,6 



ML 



M 7.7 



M 2 
Ml 



ML = V2m 2 BU (v' u /v u ) + XiAiv' s {v d /v u ) + X 5 A 5 v s {v' d /v u ) - Xfv 2 d /2 - X 2 v 2 /2 

1 



[(A1A4 + X 3 X 5 )v s v' s + (A4A5 + XiX 3 )v' d Vd] (v' u /v u ) + 
1 



M +91) + ^ 



xyj2 + xivi/2 



\{0y + to) + 



ML 



-m 2 BU + [(A1A4 + X 3 X 5 )v' s v s + (A4A5 + XiX 3 )v d v' d ]/V2 



+ V2 



;(9y + 9l) + *gl 



VuV u , 



V2m 2 BU (v u /v' u ) + X 3 A 3 v' s (v d /v' u ) + X 4 A 4 v s (v d /v' u ) 
[(A1A4 + X 3 X 5 )v s v' s + (A 4 A 5 + XiX 3 )v' d Vd] (««/<) + 2 



\{g 2 Y + gl) + 



M 2 5 = V2m 2 BD {v' d /v d ) + XiAnttvJva) + X 4 A 4 v s {v' u /v d ) - X 2 v 2 j2 - X 2 4 v 2 /2 

1 



[(A1A5 + X 3 X 4 )v s v' s + (A4A5 + XiX 3 )v' u v u ] (v d /v d ) + v d 
1 



XX/2 + \lvi/2 + 



M, 



5,6 



-m 2 BD + [(A1A5 + X 3 X 4 )v' s v s + (A4A5 + XiX 3 )v u v' u ]/V2 



+ V2 



-{g 2 Y + gi) + $gl 



VdV d , 



V2m 2 BD (v d /v' d ) + X 3 A 3 v' s {v' u /v' d ) + X 5 A 5 v s {v u /v' d ) 
[(A1A5 + X 3 X 4 )v s v' s + (A4A5 + XiX 3 )v' u v u ] (v d /v' d ) + 2 



(g 2 Y+gl) + $gl 



(v' d ? 



AL 



8,8 



V2m 2 BS {v'Jv s ) + X 4 A 4 {v' u v d /v s ) + X 5 A 5 (v u v' d /v s ) - Xjvj/2 - A§t£/2 
[(A1A5 + X 3 X 4 )v d v' d + (A1A4 + X 3 X 5 )v' u v u ] (v'Jv s ) + [25 gl] v 2 , 
X 2 v 2 d /2 + Xlv 2 u /2+[25g 2 x ] v 2 , 

M L = ~ m BS + I( A i A 5 + X 3 X 4 )v' d v d + (A1A4 + X 3 X 5 )v u v' u ]/V2 + V2 [25g 2 x ] v s v' s . 

V2m 2 BS (v s /v' s ) + XiA 1 (v u v d /v' a ) + X 3 A 3 (v' u v' d /v' s ) 

[(A1A5 + X 3 X 4 )v d v' d + (A1A4 + X 3 X 5 )v' u v u ] (v s /v' s ) + 2 [25g 2 x ] (v' s ) 2 , 
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A J 2 

M lA 


- M 2,5 


Ml,, 


= Ml, 


Ml fi 


= Ml, 




+ V2 


Mi A 


= Ml, 




+ V2 



ML 



M 
M 
M 
M. 



1.7 

2 

1,8 
2 

1.9 

2 

3,7 



Ml, 



M 2 
M 
M 
M, 



M, 



4.7 

2 

4,8 
2 

4.9 

2 

6,7 
2 

6,9 



M 2 5 = -\iA lV ' s + ZX 2 v u v d /2 + (A 4 A 5 + XiX 3 )v' u v'd + 



V u Vd, 



X\v u v d /2 + 



~\(9y +9%) + & 9l 



V u Vd, 



= -X 5 A 5 v s /V2 + V2X 2 ,v u v' d + (A4A5 + XiX 3 )v' u v d /V2 
1 



;(9y +9l) + 6g x 



v u v d , 



= -\,A,v s /V2 + V2\\v' u v d + (A4A5 + \i\ 3 )v u v' d /V2 
1 



v u v d , 



X 3 A 3 v' s + 2X 2 3 v' u v' d + (A1A3 + X,X 5 )v u v d + 2 



95) + Zgt 

Mis = -X 5 A 5 v' d + 3X 2 v u v s /2 + (A X A 4 + A 3 A 5 )« + [-10^] v u v 



(9y 



V' u v'd, 



M. 



2,7 



Xlv u v s /2+ [-IO5 2 ] v u v s , 



M| i9 = -X x A lVd /V2 + V2Xfv u v' a + (A1A4 + X 3 X 5 )v' u v s /V2 + V2 [-Wgl] v u v' s , 

Ml 8 = -X,A,v d /V2 + V2Xlv' u v s + (A1A4 + X 3 X 5 )v u v'jV2 + V2 [-10g 2 x ] v' u v s , 

-X 3 A 3 v' d + 2X 2 3 v' u v' s + (A1A4 + X 3 X 5 )v u Vs + 2 [-10 5 2 ] v' u v' s , 

M 2 a = ~X,A,v' u + 3X 2 ,v d v s /2 + (A1A5 + X 3 X,)v' d v' s + [-I5g 2 x ] v d v s , 

Ml-j = X\v d v s /2 + [-15.g 2 ] v d v s , 

M 2 g = -\iAivjV2 + V2Xlv d v' s + (A1A5 + X 3 X,)v' d v s /V2 + V2 [-15g 2 x ] v d v' s , 
M ls = -X 5 A 5 v u /V2 + V2Xjv' d v a + (A1A5 + X 3 X,)v d v' s /V2 + V2 [-lbg 2 x ] v' d v s , 
-X 3 A 3 v' u + 2Xjv d v' s + (A1A5 + X 3 X,)v d v s + 2 [-I5g 2 x ] v' d v' s . 



Neutral CP odd Higgs boson 

iPU,a 



Hi 



D 



D _ 1pD,a 



V2 ' 



V2 ' 



S a D^ (a = 1,2,3), 



/ PU.a 

C D -piM^Pj, Pi = p D ,a 

PS,a 



V2 

= 1,2, 



,9), 



(113) 



Mi 2 ! = M 2 2 = s/2m 2 BU {v' u /v u ) + XiAW s {vd/v u ) + X 5 A 5 v s (v' d /v u ) - X\v 2 d j2 - X 2 ,v 2 /2 

- [(A1A4 + X 3 X 5 )v s v' s + (A4A5 + XiX 3 )v' d Vd] (v' u /v u ), 
Ml 2 = X 2 v 2 /2 + X 2 v 2 /2 7 

M 2 3 = M 2 2 3 = -mly + [(A1A4 + X 3 X 5 )v' s v s + (A4A5 + X!X 3 )Vdv' d ]/V2, 

M 3 2 3 = v^mluK/O + X 3 A 3 v' s (v' d /v' u ) + X,A,v s (v d /v' u ) 

- [(A1A4 + X 3 X,)v s v' s + (A4A5 + X 1 X 3 )v d v d ] {v u /v' u ), 

Ml, = Ml, = V2m 2 BD (v' d /v d ) + XiAiV s (v u /v d ) + X,A,v s {v' u /v d ) - X\v 2 u /2 - X 2 4 v 2 /2 

- [(A1A5 + X 3 X,)v s v' s + (A4A5 + X 1 X 3 )v' u v u ] (v' d /v d ), 
Ml, = X 2 v 2 /2 + X 2 v 2 u /2, 

M 4 2 6 = M 2 6 = -m% D + [(A1A5 + A 3 A 4 )w> s + (A 4 A 5 + AiA 3 )u u <]/\/2, 

M 2 6 - V2m| D (V^) + ^A 3 v' s (v'Jv' d ) + X 5 A 5 v s (v u /v' d ) 

- [{XiX, + X 3 X,)v s v' s + (A4A5 + X 1 X 3 )v' u v u ] (v d /v d ), 

Ml 7 = Ml s = \^rn% s (v' s /v s ) + X,A,{v' u v d /v s ) + X 5 A 5 (v u v' d /v 3 ) - X 2 ,v 2 d /2 - X 2 ,v 2 j2 

- [(A1A5 + X 3 X,)v d v d + (A1A4 + AsAs)^^] (v'Jv s ), 
Mis = X 2 ,v 2 /2 + X 2 ,v 2 u /2, 
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Af£ 9 = Af 8 2 9 = -m2 JS + [(AiA5 + A 3 A4)w>d + (AiA4 + A3A5)w tl <]/%/2, 
M 2 9 = V2m 2 BS (v s /v' s ) + Ai^iK« rf /^) + A 3 A 3 «^M) 
- [(AiA 5 + A 3 A 4 )?;di'd + (A1A4 + X 3 X 5 )v' u v u ] (v s /v' s ), 



M l,4 


i\ r2 

= M 2,5 = 


XiA x v s + (A4A5 - XiX 3 )v u v d - X^uVd/2, 




71 f 2 




M l,6 


71 T 2 

= M 2i6 = 


A5A5U S /V2 + (AiA 3 - A 4 A 5 )?j u ?j d /v2, 


71 ,r2 

M 34 


= M 3 2 5 - 


A 4 A 4 u s /V2 + (A1A3 - A4A B )w„Ud/V2, 


M 3,6 


= A 3 A 3 i; s 


+ (A4A5 - XiX 3 )v u v d , 


?\ r 2 


71 /f2 

= M 2,8 = 


\ A f 1 ( \ \ \ \ \ 1 1 \2 /o 

X 5 A 5 v d + (A1A4 - AjjAsJv^ - X 5 v u v s /2, 


A/f 2 
M l,8 


- iW 2,7 - 




71 ,r2 

<9 


= M 2 2 9 = 


A1A1VV2+ (A3A5 - A 1 A 4 )<u s /v / 2, 




= M 3 2 8 = 


A 4 A t u d /V2+ (A3A5 - AiA 4 )w„^/V2, 


^3,9 


= A 3 A 3 u^ 


+ (A1A4 - A 3 A 5 )v„v s , 


M 4 2 7 


= Ml 8 = 


X 4 A 4 v' u + (A1A5 - X 3 XA)v' d v' s - X\vdV s /2, 


^4,8 


= ^5,7 = 


X\v d v s /2, 


^4,9 


= ^5,9 = 


XxA-lVu/S + (A3A4 - XiX 5 )v' d v s /V2, 


Ml 7 


= Ml 8 = 


X 5 A 5 v u /V2 + (X 3 X 4 - AiA 5 )w d </V2, 


Ml, 


= X 3 A 3 v' u 


+ (A1A5 - A 3 A 4 )u d w s . 



Charged Higgs boson 

H^DH+ 7 H® D H~ +3 (a = 1,2,3), 

-C D li; M*H- . H+ = ( { Ji^ )f ) = 1,2,- ■ ■ , 6), 



M'li = V2m 2 BU (v' u /v u ) + X^AWAvd/Vu) + X 5 A 5 v s (v' d /v u ) - X 2 v d - A 2 [tj 2 /2 + [v' d f] 

- [(A1A4 + X 3 X 5 )v s v' s + (A 4 A 5 + AiAs)^^^] (v' u /v u ) 

- \glK + {v' u ) 2 -v 2 d -W d ) 2 ] + \g 2 v 2 u , 

M? a = X 2 v 2 /2+^g 2 v 2 u , 

M 2 3 = M 2 3 = -m 2 BU + (A1A4 + X 3 X 5 )v s v' s /V2 + ^-j=g 2 v' u v u , 

M 2 3 = V2m 2 BU (v u /v' u ) + X 3 A 3 v' s (v' d /v' u ) + X 4 A 4 v s (v d /v' u ) - X\v 2 d - X\{v' d ) 2 

~ [(A1A4 + AsAs)^^ + (A4A5 + AiA^u^] (v u /v' u ) 

- ^ 2 2 b 2 + K) 2 -^"(^) 2 ] + ^ 2 2 «) 2 , 

M 2 4 = M 2 5 = V2m 2 BD (v' d /v d ) + XiA±v' s (v u /v d ) + X 4 A 4 v s (v' u /v d ) 

- X\vl - X\[v 2 J2 + (v' u ) 2 } - [(AiA 5 + X 3 X 4 )v s v' s + (A 4 A 5 + AiA 3 )u>„] (v' d /v d ) 

+ l9 2 2K + K) 2 -v 2 d -(v' d ) 2 } + \g 2 v 2 , 

Ml 5 = X 2 v 2 j2+\g 2 2 v 2 , 

M'l 6 = M| i6 = -m| D + (AiA 5 + A 3 A 4 )u s ^/V2 + ^=g 2 v' d v d , 

Ml 6 = V2m 2 BD (v d /v d ) + X 3 A 3 v' s {v' u /v' d ) + X 5 A 5 v s {v u /v d ) - A 2 tj 2 - A 2 «) 2 

- [(A1A5 + A 3 A4)v s ^ + (A4A5 + AiA 3 )^u„] (vd/v' d ) 
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\gl[vl + K? -v 2 -(v> d ) 2 ) + \gi{v d f, 



Ml 



Ml 



M. 



2,5 



XiAiv' s - \i[\iv u v d + X 3 v' u v' d ] + -g 2 v u v d , 



M 2A = -^glvuVd, 



Ml 6 



M 



M 



3,4 
2 

3,6 



M. 



3,5 



\ 5 A 5 v s /V2 - \ 5 [\ 4 v' u v d + X 5 v u v' d }/V2 + -^-j=glv u v' dl 
XiAiVs/V? - A 4 [A 4 ^w d + X 5 v u v' d ]/V2+ ^~^92 v 'u v <i, 



X 3 A 3 v' s - X 3 [Xiv u v d + X 3 v' u v d ] + -g\v' u v' d . 



(115) 



1-loop corrections to Higgs mass 

In order to satisfy the experimental bound for the lightest neutral CP even Higgs boson mass, the 
contributions from 1-loop corrections are important [24] [25] . We add 1-loop contributions 



AV 



1 



64tt 2 



Str 



M 2 3 



(116) 



to Higgs potentail.The dominant contributions are given by trilinear terms Y t H 3 Q 3 U 3 and kS 3 (<?i<?i + <?2<?2 + 
g 3 g 3 ). From the mass terms of squark and scalar g-quark 

-C D [m 2 Q3 + (Y t (HV) ) 2 ] \U 3 \ 2 + \m 2 + {Y t {H^f) 2 \ \U C 3 \ 



3 T 

/ \21 /| „c|2 



+ [m 2 g + (kv' s ) 2 ] (\ 9l \ 2 + \g 2 \ 2 + |. 93 | 2 ) + K + (K) 2 ] (\9 c i\ 2 + \9 C 2 \ 2 + l<? 3 T) 
+ [Y t A t (HV)° + Y t X 3 v> s ((H°fy] U 3 U C 3 
+ [kA k v' s + kX 3 ((HU) (H°)°y] { 9l gl + g 2 g c 2 + g 3 g c 3 ) + h.c. 



m 2 Q3 + (Y t v' u ) 2 Y t A t v' u + Y t X 3 v' s v' d 



' \2 



Y t A t v' u + Y t X 3 v' s v' d m 2 +{Y t v'J 



(kv' s 



kA k v' s + kX 3 v' u v' d 



kA k v' s + kX 3 v' u v' d m 2 + (kv' s ) 



I \2 



u 3 

my 

gi 

(<??)* 



(117) 



M 
M 



T,± 
2 



mass eigenvalues are given by 

I 
2 
1 

L g,± = 2 
For simplicity, we assume 



m Q 3 + m k + m<) 2 ± y/(m 2 Q3 ~ ™ 2 a S ) 2 + ±Yt 2 (AtV u + X 3 v' s v' d ) 2 
m 2 + ml + 2(kv' s ) 2 ± yj (m 2 - m 2 .) 2 + 4k 2 (A k v' s + X 3 v' u v' d ) 



m 2 g = ml 



16TeV^ 



A t = A k 



O.OTeV, 



then Eq.(118) is rewritten by 



M. 



M 



T,± 
2 



9,± 



m Q + ( Y tv' u ) 2 ± Y t X 3 v' s v' d , 
m Q + (H) 2 ± kX 3 v' u v'd- 



As potential minimum condition is modified as 

d(V + AV) 



and we must add the terms 



dX 



AM 2 , 3 



AA/ 6 2 6 



AM 



.3.6 



= 0, X = (v' u ,v' d ), 

1 d 2 AV 1 dAV 
2d(<? '" 2«) d(v' u Y 
1 d 2 AV 1 dAV 
2d{^P~W^ d(v d y 
1 d 2 AV 
2d(v> u )d(v' d y 



(118) 
(119) 

(120) 
(121) 



(122) 



17 



to the neutral CP-even Higgs boson mass. We fix renormalization point as 

A = 4TeV. 



Chargino 

C 



3 Xi M ijXj +h.c, 

.v _ ( K' 
hi 



M = 



Neutralino 



u la 
+ 



X, 



h°, 
"d 



(a = 1,2,3), 



—iw 



^2 ~F 
V2 



( \iv' s X 5 v s /V2 g 2 v u /V2 \_ 

X t v' s X 5 v s /V2 g 2 v u /V2 

XiVs/y/2 A 4 u s /v / 2 A 3 ^ g 2 v' u 



V g 2 v d /V2 g 2 Vd/V2 gW d 



Mo 



J 



h.c. 



X 



M 



( (K)a \ 
(h°d)a 
Sa 

v * J 

( M ud M, 

M du M ds M dx 

M su M sd M sX 

\ M Xu Mxd M Xs M X \ ) 




(123) 



(i,j = 1,2,3,4), 



(a =1,2, 3; i,j 



1,2,. 



,12), 



(124) 



ltd 



M„ 



A/, 



M, 



sA 



A/aa = 



AC, 



A 5 w s /x/2 \ 
X 5 v s /V2 
Mv' s J 

X lVd /V2 \ 
Xiv d /V2 = Mj u , 
hv' d J 

Xiv u /V2 \ 
X lVu /V2 



X lV ' s 
Aii^ 
X 4 v s /V2 X 4 v s /V2 

X 5 v' d 
X 5 v' d 
A 4 Ud/\/2 XiV d /V2 

XW U 
X 4 v' u 
X 5 v u /V2 X 5 v u /V2 X 3 v' u J 

gYV u /2 -g 2 v u /2 -2g x v u 
g Y v u /2 -g 2 v u /2 -2g x v u 
g Y v' u /V2 -g 2 v' u /V2 -2V2g x v' u 

-g Y v d /2 g 2 v d /2 -3g x v d 
-g Y v d /2 g 2 v d /2 -3g x v d 
-g Y v' d /y/2 g 2 v' d /V2 -3V2g x v' d 

5g x v s 
5g x v s 
5V2g x v' s 

-M Y 
-M 2 
-M x 



M, 



Ml, 



= M 



Xu i 



(125) 



B Mixing matrices and mass eigenvalues 

Mass eigenvalues (m^ : GeV) and diagnalization matrix U = (u\, u 2l ■ ■ •) are given as follows. 
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Neutral CP-even Higgs 

For Higgs bosons, the diagonalization matrices are defined as 

(M 2 )' = U T M 2 U = diag(m 2 , m 2 , ■ ■ ■). (126) 



i 


1 


2 


3 


4 


5 


6 


7 


8 


9 


Ui 


-0.698 


-0.113 


-0.0006 


-0.613 


0.329 


-0.115 


0.0006 


0.003 


-0.059 




0.698 


0.113 


0.0006 


-0.613 


0.329 


-0.115 


0.0006 


0.003 


-0.059 













-0.384 


-0.808 


-0.036 


0.029 


-0.009 


-0.445 




0.0015 


-0.005 


-0.707 


-0.0012 


-0.008 


0.006 


-0.706 


-0.0008 


-0.030 




-0.0015 


0.005 


0.707 


-0.0012 


-0.008 


0.006 


-0.706 


-0.0008 


-0.030 













-0.273 


-0.360 


-0.031 


-0.033 


0.054 


0.889 




0.113 


-0.698 


-0.0056 


0.118 


-0.025 


-0.693 


0.006 


0.071 


-0.002 




-0.113 


0.698 


0.0056 


0.118 


-0.025 


-0.693 


0.006 


0.071 


-0.002 













-0.002 


0.014 


0.101 


0.00009 


0.993 


-0.052 


mi 


91.50 


537 


2005 


121.96 


152.96 


539 


2007 


1018 


1425 



(127) 



Neutral CP-odd Higgs 



i 


1 


2 


3 


4 


5 


6 


7 


8 


9 




-0.703 


0.073 


-0.040 


0.004 


0.0005 


0.700 


0.074 


-0.002 


-0.059 




0.703 


-0.073 


-0.040 


0.004 


-0.0005 


0.700 


0.074 


-0.002 


-0.059 










-0.889 


0.086 





-0.088 


-0.009 


-0.029 


-0.044 




0.0014 


0.008 


0.008 


-0.00008 


0.707 


0.0003 


0.008 


-0.706 


0.030 




-0.0014 


0.008 


0.008 


-0.00008 


-0.707 


0.0003 


0.008 


-0.706 


0.030 










0.448 


-0.004 





-0.049 


-0.005 


-0.033 


-0.892 




-0.073 


-0.703 


-0.006 


-0.070 


0.008 


0.074 


-0.700 


-0.008 


-0.002 




0.073 


0.703 


-0.006 


-0.070 


-0.008 


0.074 


-0.700 


-0.008 


-0.002 










-0.079 


-0.991 





-0.012 


0.099 


-0.001 


-0.034 




112.12 


532 








2005 


145.53 


535 


2008 


1423 



(128) 



Charged Higgs 



i 


1 


2 


3 


4 


5 


6 


Ui 


0.707 


-0.000055 


-0.703 


0.041 


-0.001 


0.059 




-0.707 


0.000055 


-0.703 


0.041 


-0.001 


0.059 










0.089 


0.893 


-0.029 


0.441 




0.000055 


0.707 


-0.0018 


-0.008 


-0.706 


-0.030 




-0.000055 


-0.707 


-0.0018 


-0.008 


-0.706 


-0.030 










0.049 


-0.447 


-0.032 


0.893 


rtii 


90.70 


2005 


130.02 





2007 


1422 



(129) 



Neutralino 

For neutralinos, the diagonalization matrix is defined as 
U T MU = diag(mi,m2, ■ ■ ■ Jnia). 



19 



i 


1 


2 


3 


4 


5 


6 


Ui 


0.234 


-0.496 


0.446 


-0.248 


-0.475 


0.067 




-0.234 


0.496 


-0.446 


-0.248 


-0.475 


0.067 













0.064 


0.148 


0.031 




0.209 


0.504 


0.450 


-0.217 


0.497 


-0.055 




-0.209 


-0.504 


-0.450 


-0.217 


0.497 


-0.055 













0.028 


-0.105 


-0.066 




— U.0O4 


— u.ui / 


n oiq 
U.olo 


n £91 


ft H9Q 
— U.Uzo 


— u.uuz 




0.634 


0.017 


—0.313 


0.621 


—0.023 


—0.002 




n 
u 


n 
u 


n 

u 


u.uou 


u.uuo 


n nnnn7 




n 
u 


n 
u 


n 


n nno 
u.uuy 


U.U04 


u.yoo 




n 
u 


n 
u 


n 


n m a 

U.U14 


n 1 1 o 


n 991 




n 
u 


n 
u 


n 


U.UUU4 


n nm 

U.UUl 


n fin'? 


m. 


A 1 09 


— loU 


1 79 


/i/i 

— 44.00 


— in 


1 77 
— Ill 


i 


7 


8 


9 


10 


11 


12 


Ui 


0.062 


-0.430 


-0.099 


-0.099 


-0.020 


0.002 




0.062 


-0.430 


-0.099 


-0.099 


-0.020 


0.002 




0.039 


0.132 


-0.689 


-0.682 


-0.121 


-0.003 




-0.049 


-0.439 


0.062 


-0.063 


-0.013 


-0.003 




-0.049 


-0.439 


0.062 


-0.063 


-0.013 


-0.003 




-0.077 


0.085 


0.692 


-0.686 


-0.138 


-0.050 




— 0.002 


a o o a 
— 0.330 


a aai 
0.001 


a aao 
— 0.003 


—0.053 


n A A 7 

0.04/ 




-0.002 


-0.330 


0.001 


-0.003 


-0.053 


0.047 




-0.0001 


0.051 


0.014 


0.114 


-0.728 


0.669 




-0.241 


0.007 


0.063 


-0.014 


-0.002 


-0.001 




-0.960 


-0.012 


-0.118 


0.024 


0.003 


0.002 




0.003 


-0.0009 


-0.032 


0.153 


-0.655 


-0.739 


mi 


-194 


157 


832 


820 


945 


-1144 



Chargino 

For the charginos, the diagonalization matrices arc defined as 



X 

uTmu. 



U 4 



u-(x~y> x + = u + (x + y, 

= diag(mi, rri2, Jn.3, m^) = (130, 111, 193, 830), 

/ 0.707 0.685 0.104 0.140 \ 

0.707 0.685 0.104 0.140 

-0.209 0.066 0.976 

\ 0.131 -0.987 0.095 / 



(130) 



U- 



( -0.707 -0.696 

0.707 -0.696 

0.140 

\ -0.102 



0.085 0.088 \ 

0.085 0.088 

0.128 0.982 

-0.984 0.143 J 



(131) 



These mass eigenvalues are consistent with the experimental mass bounds [18] 



Charged Higgs 
Neutral CP-even Higgs 
Neutral CP-odd Higgs 
Chargino 
Neutralino 



m > 79.3GeV, 
m > 114.4GcV, 
m > 93.4GcV, 
m > 94GeV, 
m > 46GeV. 



Note that p' 3 , p' 4 , (Hf)' are Nambu-Goldstone boson which are eaten by gauge bosons. 



C The lifetimes of exotic quarks 

As the R-parities of exotic quarks are odd, at least there must be one sfcrmion which is lighter than exotic 
quarks, to make them unstable. Now we assume the right handed slepton Ef is lighter than 2TeV and the 
other sfcrmions are heavier than ATeV. For simplicity, we assume there is no mixing between El and Li. 
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Through the non-renormalizable interaction 

C D - E ^i"^' 1 + + ^,3)^< + h.c, (132) 

the exotic quarks ipg,i~3 can decay into u\ and E\ , where Cj are O(l) coefficients. The lifetimes are estimated 
as follows 

from which we must put Ci ~ 4 in order to satisfy the cosmological constraint for exotic particle, t < O.lsec. 
The interaction Eq.(132) comes from 

W D E jj^sC^i + $2.92 + *§sa)£?E/?, (134) 
which may contribute to LF decay through 

C D ~ E Jp a j ( t ) LF^ g ,jE c 1 u c l + h.c, (135) 
i,0 P 

where otj arc given by linear combinations of the flavon mixing parameters. In this paper, we assume ctj are 
small enough and this interaction does not give sizable contribution to LF decay width. 
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